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ABSTRACT 

Spectroscopic galaxy-galaxy lens searches are presently the most prolific method of identifying strong 
lens systems in large data sets. We study the probabilities associated with these lens searches, namely 
the probability of identifying a candidate with rogue [O II] emission lines in a galaxy's spectrum, 
and the probability that the candidate will show features of strong lensing in follow-up photometric 
observations. We include se lection effects uniqu e to spectroscopic data, and apply them to the Sloan 
Lens ACS (SLACS) survey {Bo lton et alj|206"q ). The most significant selection effect is the finite size 
of the spectroscopic fiber which selects against large separation lenses and results in a non-monotonic 
dependence of the rogue line probability on velocity dispersion. For example, with the 3 arcsec 
diameter SDSS fiber and 2 arcsec FWHM seeing, we find that the probability that a given LRG has 
a rogue [O II] line in its spectrum decreases with velocity dispersion from 150 km/s to 300 km/s and 
then increases up to 400 km/s for a given source size. The total probability for observing a rogue 
line in a single survey spectrum is ~0. 9-3.0%, and the total lensing rate is ~0.5-1.3%. The range is 
due to uncertainties in the physical size of [O II] emission regions, and in the evolution of the [O II] 
luminosity function. Our estimates are a factor of ~5 higher than the results of the SLACS survey, 
a discrepancy which we attribute to the SLACS requirement that multiple rogue lines be observed 
simultaneously. 

Subject headings: gravitational lensing - surveys - galaxies: statistics 



1. INTRODUCTION 

Spectroscopic gravitational lens searches have begun 
to yield a remark able number of strong galaxy-galaxy 
(g-g) lens systems (|Bolton et al.ll2004 120061; IWillis et all 
2006). These finite source lenses promise both new phys- 
ical insights and new phenomenology. The extended 
images provide extensive constraints on the lens poten- 
tial, especially on the radial density profile, which is still 
the main systematic uncerta inty in lensing constraints 
on the Hubble c onstant (e.g. jKeeton fc K ochanckl ll997t 
Kochanck 2002). In present surveys, limits on the source 
rcdshift range mean the lenses that are found typically 
have images that a ppear well within the e ffective radius 
of the lens galaxy (Koopmans et al. 2006). This makes 
g-g lenses ideal for probing the inner regions of distant 
elliptical galaxies. In addition, current selection effects 
favor star-forming source galaxies, which opens the excit- 
ing possibility of observing multiply-imaged sup ernovae 
(|Oguri fc Kawanoll2003t iDobler fe Keetonil2006| i. 

The basic premise behind spectroscopic lens searches 
is to mine large samples of galaxy spectra looking for 
"rogue" emission lines that originate fro m background 
galax i es at small impact parameter (j Warren et alJ 
119961: IWillial |2000D . This tech nique is complemen- 
tary to photometric s e arche s (e.g.. ICabanac et all 120071 : 
iKubo fc DeirAntoniol 120071 ) which look for strongly 
lensed, arc-like features in imaging data. 
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Among several recent spectroscop ic searches 

(|Bolton et all l200i I200& IWillis et alJ l2006h . the 
most prolific has been the S loan Lens A C S (SL ACS) 
survey. 6 F o r this survey, iBolton et alJ (|2004h and 
IBolton et ail (|2006l ) mined a catalog of 50,996 Sloan 
Digital S ky Survey (SPSS) L uminous Red Galaxy 
(LRG, see Eis enstein et al.ll2001l ) spectra for rogue [O II] 
3727 emission lines, and found ~50 candidates. The 
addition of later SDSS data releases as well as spectra 
from the MAIN galaxy sample increased the number of 
candidates to ^200, with a similar l-in-1000 incidence. 
Follow-up observations have subsequently confirmed 
70-80 new g-g lenses from among these candidates. 
SLACS data have been u sed to place the lens galaxies on 
the f undamental plane (|Treu et all 120061 : IBolton et al.l 
2007), to constrain the redshift evolution of the densit y 
profiles of elliptical galaxies (jKoopmans et all l2006h . 
and to trace the den sity profiles out to very large radii 
(|Gavazzi et alll2007f) . 

Given that galaxy-galaxy lenses are already numerous, 
and will become incr e asingly common in lar ge surveys 
(|Marshall et all 120051 : iMoustakas et all I2007T) . a sound 
statistical analysis of the expected incidence of g-g strong 
lensing is warranted. In trad itional analyses of t he statis- 
tics of lensed quasars (e.g.. iTurner et al1l!984h the pri- 
mary statistical question is, "what is the probability that 
a given source is lensed?" By contrast, in g-g lens statis- 
tics the question is different, viz. "what is the probability 
that a given galaxy is a lens?" In this paper we formulate 
a general statistical analysis applicable to spectroscopic 
g-g lens searches, and apply it to the SLACS sample to 
estimate the total number of rogue emission lines in the 
survey, and the actual number of strong lens systems that 
should be confirmed by follow-up observations. Our re- 

6 www.slacs.org 
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suits for SLACS will help assess the completeness of that 
survey, while our general conclusions will (we hope) be 
useful in the design of future spectroscopic lens searches. 

Except where noted, throughout this paper we assume 
a flat cosmology with VI m = 0.3, Vl\ = 0.7, and Hq = 70 
km/s/Mpc. 

2. PROBABILITY FOR GALAXY-GALAXY LENSING 

Figure [1] shows a schematic representation of the lens- 
ing geometry. Let the galaxy be described by a singular 
isothermal ellipsoidal (SIE) mass distribution with pa- 
rameters G — (zi, a, e, 7, 7 ) where zi is the redshift, a 
is the velocity dispersion, e is the ellipticity, 7 is the ex- 
ternal shear, and is the angle between the ellipticity 

and shear. The probability Pq(G) that this galaxy is a 
lens is equivalent to the probability that there is a source 
within the "Einstein cone" of the galaxy. Here we define 
the Einstein cone to be the region behind the galaxy in 
which a source is strongly lensed. For point sources, the 
Einstein cone is the same as the multiply-imaged region, 
but we will refine the definition shortly to incorporate 
complexities from extended sources. Note that the cross 
section of this region need not be circular, and its size 
does not grow linearly with distance, so the volume is 
not strictly conical; but we believe the terminology is 
convenient and attractive. 

If the number density of sources brighter than flux S 
as a function of redshift is given by n s (z s ,S), then the 
lensing probability is 



P G (G)= n s (z s ,S)dV 

roo d y 



(1) 



dz«dVL 



dz. 



n s (z Sl S) du, 



where V e i n is the volume of the Einstein cone, zi and 
z s are the lens and source redshifts respectively, and 
= Vt(z s ) is the solid angle subtended by the cone 



VI 



at z s . The integral over VI is actually an integral over 
u, the angular coordinates in the source plane, and for 
now we consider the source plane integral to extend over 
the multiply- imaged region. Finally, it is natural to do 
the integral in comoving coordinates, but to express dis- 
tances as angular diameter distances. So we write the 
comoving volume element as 



dV 
dz*dVl 



c_ (l + z s ) 2 Dl 
H 



E(z s 



(2) 



where E(z) — [Vl M (\ + z) 3 + SIa] 1 / 2 , and write factors 
of (1 + z s ) explicitly allow us to keep D s as the angular 
diameter distance to the source. 

In practice, we are interested in the number density 
of sources whose observed flux is above a survey's flux 
limit Sq. This implies Si x fi > So where Si is the 
source's intrinsic flux and /x is the lensing magnification. 
Therefore, the relevant number density for eq. ([1]) is 

/>oo 

n s (z s ,S) =n s (z s , So/fi) = / $(L,z s )dL, (3) 

J L /p, 

where <fr(L,z s ) is the source luminosity function at red- 
shift z s . Here Lo = 4ir(l + z s ) A D 2 s So is the luminosity 
corresponding to the flux limit So, and factors of (1 + z s ) 



again appear so that we may keep D s as an angular di- 
ameter distance. The fact that the lower limit of integra- 
tion depends on \i means that the integral automatically 
incorporates lensing magnification bias. 

With point sources the definition of a lens relies on 
image multiplicity: any source with multiple images is 
said to be strongly lensed. With extended sources the 
situation is more complicated. A source lying just out- 
side the caustics might be distorted enough to be labeled 
a lens even if there is just one image. A source lying 
astride a cusp or fold caustic may exhibit a single arc 
comprising two or three merged images, with counter- 
images that may or may not be bright enough to be 
detectable. The point is that identifying an object as 
a lens may depend on some qualitative interpretation of 
the morphology. Since the interpretation depends on dis- 
tortions of the image(s), which are related to the lensing 
magnification, we attempt to quantify the labeling of ex- 
tended lenses through a magnification cut. Specifically, 
we label an object a lens if the lensing magnification \i 
exceeds some threshold set by \i cu t- (Our choice of [i cut is 
discussed below.) We then take the source plane integral 
in eq. fl} to extend over the region in which /x > fx C ut- 

The lensing magnification depends not only on the 
source position u, but also on the parameters of the lens 
galaxy. Therefore eq. ([1]) represents the lensing probabil- 
ity for a particular galaxy. To describe a population of 
galaxies, in principle we want to average over some ap- 
propriate distributions of zi, a, e, 7, and 7 . In practice, 
if we have a set of N observed galaxies, each of which is 
described by the parameters Gi, we can do the average 
explicitly: 



1 N 

P G = ^P G (G i ). 



(4) 



i=l 



In the following sections we apply this Pq calculation to 
the SLACS sample, specifically incorporating parameters 
and selection effects appropriate to that survey. 

3. SLACS SURVEY PARAMETERS 

3.1. Spectroscopic Selection Effects 

When discussing SDSS spectra, it is crucial to account 
for the finite size of the spectroscopic fiber. The strength 
of a rogue emission line is directly related to how much 
flux from the source galaxy falls within the fiber. If a 
lens galaxy is massive, the lensed images may be pushed 
outside of the fiber, so the rogue line may be weak or 
even absent despite the fact that the system is a lens. 
Figure [3 shows a sample image configuration both with 
and without this "fiber cut" taken into account. (This 
figure ignores the effects of seeing, which are discussed in 

El) 

The diameter of the SDSS fiber is 3.0 arcsec, while an 
SIE galaxy lens with redshift zi = 0.2 and velocity disper- 
sion a = 250 km/s (roughly the mean LRG redshift and 
velocity dispersion from the SLACS survey) has an Ein- 
stein radius of ~1.0-1.2 arcsec depending on z s . Thus, 
the fiber cut may create a significant bias against large 
separation lenses in the SLACS sample. There is also a 
bias against large sources (for a given source flux). Of 
course, with follow up observations the full image con- 
figuration will be observed, but it is the fiber flux that 



3 



Observer 



Galaxy 



Multiply 

imaged 

region 



thickness = dl 

Fig. 1. — Lensing geometry for spectroscopically-selected galaxy-galaxy lenses. The probability that the foreground galaxy is a lens is 
the probability that there is a source within the galaxy's "Einstein cone." (The multiply-image region behind the galaxy is not strictly a 
cone, because the cross section need not be circular and does not grow linearly with distance; but the terminology is attractive.) 
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Fig. 2. — Selection effects due to the finite SDSS spectroscopic fiber size. (Left) The solid curves show a sample image configu- 
ration for a lens system with lens parameters (z;, cr, e, 7, </> 7 ) = (0.2, 300 km/s, 0.3, 0.05, 0.0) and source parameters (z s , R BTC , no , v ) = 
(0.6,0.88 kpc, 0.15", 0.15"). The dotted circle denotes the SDSS fiber. (Right) The solid curves show only the portion of the images that 
fall inside the fiber. The fiber cut reduces the integrated flux within the fiber, creating a selection bias against large separation lenses and 
large source sizes. (This figure ignores the effects of seeing, which are discussed in £14.21 ) 



determines whether a system is identified as a lens can- 
didate in the first place. 

The finite wavelength range of SDSS spectra places an 
upper limit on detectable source redshifts. The qual- 
ity of spectral noise modeling is another important fac- 
tor, since imperfect sky subtraction can leave emission 
line residuals that are not modeled by the eigenspectr a 
used to fit LRGs ( see Figure 1 of iBolton et al.ll2004f ). 
iBolton et al.l (|2004f ) were careful to account for imper- 
fect sky subtraction at long wavelengths, allowing them 
to probe deep into the 7000-9000 A range. A third 
factor is the SLACS selection criteria: LRG spectra 
were required to exhibit not only blended [O II] 3727 
lines, but also two longer- wavelength "secondary" fea- 
tures. The longest-wavelength secondary feature was 



[O III] 5007, while the shortest was H/3 with rest wave- 
length 4863 A. These lead to an upper limit on the 
source redshift of z s , max - 9200A/5007A - 1 = 0.84 
and z s , max ~ 9200A/4863A - 1 = 0.89, respectively. In 
our calculations we therefore take z S!max = 0.9. Finally, 
to e nsure a significant l ensing probability in their sam- 
ple, IBolton et al.l (|2004h only searched for sources with 
velocities more than 5000 km/s behind the LRG, cor- 
responding to a lower limit on the source redshift of 
Zs,min = zi + 0.017. These constitute the limits of the 
z s integration in eq. {TJ) . 

The last i mportant spec t ral par ameter is the flux limit. 
Figure 1 of IBolton et alj (|2004h shows the typical la 
noise spectrum. For wavelengths <7200 A the wave- 
length dependence is small, so for simplicity we take 
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TABLE 1 
Source O II LF Parameters 





log I/* 


a 


P 


(10 -3 Mpc" 3 ) 


(erg/a) 






4.09 


42.34 


-1.15 





2.16 


42.31 


-1.13 


1 


0.90 


42.37 


-1.17 


2 


0.62 


42.18 


-1.09 


3 



Note. — Best-fit Schechter func- 
tion parameters for the O II lumonisity 
function in the redshift range 0.3 < 
z < 1.5. We h ave c onverted the data 
of IHogg et all [fl998h from OCDM to 
ACDM, and we have considered differ- 
ent possibilities for the number evolu- 
tion parameter f3 (see text). 

the noise floor to be constant in wavelength, which 
also means consta nt in source redshift. Including the 
iBolton et al.l (|2004D requirement that secondary emission 
features have signal-to- noise ratios greater than 3.0, we 
take So = 6.0 x 10~ 17 ergs/s/cm 2 . 

3.2. Source Population: II Luminosity Function 

From eq. ([3]) it is clear that we must specify the lumi- 
nosity function (LF) of sources, in order to account for 
the flux limit and the magnification bi as. For SLACS, the 
primary selection is on [O II] line flux. IHogg et all (fl998h 
give the [O II] LF for the redshift range 0.3 < z < 1.5 
(see their Fig. 7), which covers the range of source red- 
shifts accessible in the SLACS sur vey. Howeve r, we make 
two refinements to the LF. First, IHogg et alJ quoted the 
LF for an OCDM cosmology, but we prefer to work in 
ACDM; luminosities and volumes both need to be ad- 
justed. Second, since [O II] emission is thought to trace 
the star form ation rate, the LF may vary substantiall y 
with redshift (|Kennicuttlli992t iGlazebrook et al.ll2003D . 
We include the possibility of number evolution by mod- 
eling the LF as an evolving Schechter function, 



dL = n*{l + z s f 



-L/L, 



dL 



(5) 



We make both adjustments using the following tech- 
nique. We first choose a value for the evolution param- 
eter p. We postulate a set of Schechter function param- 
eters (n*,i*,a) to specify the ACDM LF, <&acdm- We 
draw from this LF to generate a mock sample of sources 
in an ACDM universe. We then imagine "observing" 
these sources, interpreting them using an OCDM cos- 
mology, and deriving the effective OCDM LF, ^qcdm- 
We compare ^nn rw with the OCDM LF presented by 
IHogg et all (| 19981 ) to see how well they match. We then 
repeat this process for many values of (n* , L* , a) and 
choose the valu es that prov ide the best match between 
$0 CDM and the lHogg et aLl data. This gives us the best- 
fit ACDM LF, for our particular choice of the evolution 
parameter (3. Finally, we repeat the entire analysis for 
different values of [3. The resulting LFs are summarized 
in Table [TJ 

With an LF of the form eq. (J5]) , the luminosity integral 
in eq. ([3]) can be evaluated, 



/>oo 

/ *{L,z.) dL = n* (1 + z s f T 

J Lo/fJ. 



where T is the incomplete gamma function. Recall that 
for a flux-limited survey, Lq = 4tt(1 + z s ) 4 D^Sq depends 
on redshift. 

3.3. Lens Population: SDSS LRG Sample 

The initial sample analyzed bv lBolton et al.l (|2004h in- 
cluded ~51,000 LRG spectra obtained by SDSS between 
2000 March 5 and 2003 May 27. To obtain proper statis- 
tics, we must include appropriate distributions of zi, a, 
e, 7, and 4> 7 for this sample (see eq. EJ. The velocity dis- 
persion function dn/da of the full SDSS el l iptica l galaxy 
catalog has been measured bv lSheth et alJ (|2003l ). Their 
analysis corrected for various selection effects in order 
to recover the intrinsic distribution dn/da. However, we 
wish to include the selection effects since our goal is to 
estimate how many rogue emission lines should have been 
found in the actual SDSS data. 

To do this, and also to account for d istributions of zi 
and e , we randomly select 800 LRGs fsee lEisenstein et alJ 
120011 ) observed between 2000 March 5 and 2003 May 27 
and flagged as GALAXY-RED by the SDSS photomet- 
ric pipeline (|Lupton et al.ll200"l[) . 7 The number of LRGs 
was chosen to be computationally tractable, and we have 
verified that it is a sufficiently large sample to yield ac- 
curate statistics (see ^5.5|) . Choosing from the sample 
of observed LRGs automatically incorporates all of the 
same selection effe cts as the sample fro m which SLACS 
was drawn. As in Bol ton et al.l (|2004f ). we restrict the 
LRG redshift range to 0.15 < zi < 0.65. 8 We use r-band 
ellipticities from de Vaucouleurs fits in the SDSS photo- 
metric pipeline (jLupton et al.ll200lT ). These ellipticities 
describe the light while what we really need is the mass, 
but there is evidence that the mass and light elliptici- 
ties follow similar distributions (|Rusin fc Tegmarkl[2001t 
iHevl et alJH99llNaab fc BurkerlSooajT " 

We assign each galaxy a random shear amplitude 
drawn from a lognormal distributi on centered on 7 = 
0.05 with dispersion 0.2 dex (see [Holder fc Schechter] 
l2003f ). and a random shear angle 4> 1 € [0,27r]. 

4. METHODS 

Our formula for the total lensing probability explic- 
itly includes two integrals over the source redshift z s and 
the source position u (see eq.[T]). There is a third integral 
that enters implicitly: an integral over the image plane to 
compute the magnification of an extended source. The 
integral over source redshift is straightforward to com- 
pute numerically, but the 2-D integrals over the image 
and source planes require more care. 

4.1. Semi-Analytic Image Plane Integration 

To calculate the magnification /i for a given source 
and lens, we extend th e analytic method developed in 
iDobler fc Keetonl (|2005l ) to include finite source lensing 
by isothermal ellipsoids (SIEs) in an external shear field. 
The SIE density profile has been used quite successfully 



7 For detail s related to SPSS, see lYork et al.l J200(J for a techni- 
cal summary , Gunn ct al. (1998) fo r issu es rel ated to the camera , 



1 



a 



Lq 
L^\i 



, (6) 



IFukugita et all il998\ . IHogg et all J2001T) . and ISmith et all pOOl) 

for a discussion of the photometric system and calibration, and 
IPier et all l|2003f ) for details rel ated to astrome t ric ca libration. The 
tiling procedure is described in lBlanton et al. (2003). 

8 The distributions of galaxies in redshift and velocity dispersion 
are shown in Figure [9] 
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to model not only quasar le nses but also the extend ed im- 
ages seen m SLACS lenses (|Koopmans et al.ll2006D . The 
lens equation is 



_ ( 1 — 7 cos 2<\> 1 —7 sin 20 7 
—7 sin 2<j)~ 1 I + 7 cos 20 7 



x-a(f), (7) 



where x = (r cos 9, r sin 9) are image plane coordinates. 
The two co mponents of the deflection angle for an 
SIE l ens are (jKormann et al.lll994l : iKeeton fc Kochaneld 
1998) 



b'q 



b'q 

7^7 



where q = 1 



Q{0) = 
b' = 



tan -1 [Q(9) cos( 



tanlT 1 [Q(0) smt 



1/2 



q 2 cos 2 1 
bn 



sin 



2K(1 + q- 2 ) ' 



(8) 
(9) 

(10) 
(11) 



and K is the complete e lliptic integral of the first kind 
(see iHuterer et alll2005l ). The Einstein radius b of the 
galaxy is related to its velocity dispersion by 



■ c) D s 



(12) 



where D s and Di s are angular diameter distances from 
the observer to the source and from the lens to the source, 
respectively. 

The idea behind our analytic method is to parame- 
terize the source boundary by a circle: (u, v) — (uq + 
i? src cos A, v + i? src sinA) for A e [0, 2%], Plugging this 
into the lens equation yields 



uq + R S1C cos A = r T_ — a a 
Vq + R SIC sin A = r T + 



(13) 



with 



T ± = (1 ± 7cos20 7 ) cos6» - 7sin2(/> 7 sin0. (14) 



We square and add the two equations in eq. (|13[) to elim- 
inate A. Since a Xl a y , and r± are all independent of 
r, we obtain a simple quadratic equation for r that we 
can solve to find the following analytic expression for the 
image boundary as a function of 9: 



r±(9) = 



B ± VB 2 - AC 
A 



(15) 



where 



A = Y\ +T 2 _, 



B = T-(a x + u ) + T + (a y + v ), 
C=(a x + u ) 2 + (a y +v ) 2 - R^. (16) 
The total magnification of an extended source is then 
Total image area 1 



A*tot 



Total source area 



27ri? 2 rc j M 



(rl-rl 



(19 
(17) 

where M. is the region of 9 over which r±{9) is real and 
positive. We can also impose the SDSS fiber cut very 



F--S.17 



L = 5.50 




m 



tm - 0.07 





Fig. 3. — Sample image configurations for on-axis (top) and 
off-axis (bottom) lensing geometries both with (right) and without 
(left) seeing. Here i? c i n = 0.9 arcsec and R BTC = 0.07 arcsec, the 
source flux is normalized to unity, the circle represents the = 
1.5 arcsec SDSS fiber, and the seeing is 2 arcsec (FWHM). For on- 
axis sourcees, the fiber flux is preferentially decreased due to seeing, 
while for off-axis sources the fiber flux is preferentially increased. 



simply as follows. Let r^b = 1.5 arcsec be the fiber ra- 
dius. If we define the "fiber magnification" to be the 
total flux within the fiber divided by the total flux of the 
source, we can compute this as 



1 



Mfib 



2ttR 2 



M 



[min(r + , r fib ) 2 - min(r_ , r fib ) 2 ] d9 

(18) 

We emphasize that in both cases the solution for the im- 
age boundary is completely analytic, which allows us to 
reduce the 2-D image plane integration to a 1-D integral 
along the image boundary. 

It is useful to note explicitly where the various system 
parameters enter into the magnification calculation. The 
lens and source redshifts and the galaxy velocity disper- 
sion all enter through the Einstein radius b. The shear 
parameters appear in T±, while the ellipticity is buried 
in the deflection components a x and a y . 

4.2. Seeing Effects 

The above treatment of the image plane integration ne- 
glects the effects of atmospheric seeing. For the SLACS 
sample, the seeing is typically ~2 arcsec (FWHM), which 
is an appreciable fraction of the SDSS fiber diameter. 
Seeing can either add flux to the fiber from images out- 
side, or remove flux from the fiber from images that are 
inside. Which of these effects dominates depends on the 
configuration of images, as shown in Figure [3l 

We handle the effects of seeing in the following way. 
Let J(x) be the surface brightness distribution of the 
lensed image (in the absence of seeing). Then the 
smeared surface brightness distribution is 



J'(x) = J G(x-x') J(x') dad 



(19) 
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where G represents the PSF, which we take to be a Gaus- 
sian with FWHM 2 arcsec. Specifically, G(x — x') is the 
flux at x when the PSF is centered at x'. The fiber flux 
is then 

F Rh = [ J'(x) dx 

JRh 

= [ dx [ dx G(x - x') I(x) 



[ dx' I(x') f dx G(x - x 1 



(20) 



The term in square brackets is the fiber flux of a Gaussian 
centered at x 1 . Since the Gaussian PSF and the fiber are 
both circular, this term can only depend on the distance 
of the center of the Gaussian from the origin, r = \x'\. 
Hence we write this factor as Gf^(r), an d then rewrite 
cq. (HOI) as 



Fab = 



dr i 



2,T 



dS I(r,6) Gfib(r). 



(21) 



If we take the source flux to be unity, the integral ac- 
tually gives the fiber magnification /xgb- In this case, 
the source surface brightness is l/nR 2 rc , and since lens- 
ing conserves surface brightness we have I(x) = 
within the image boundaries, and outside. Since the 
image boundaries are given by r± (8) from eq. (| 15(1 , in the 
end we can write the fiber magnification in the presence 
of seeing as 



1 



dd 



r+(e) 



-(6) 



dr r Gfib(r), (22) 



which is a generalized version of eq. ([Tall . The "fiber 
Gaussian" factor Gab(r) must be computed numerically, 
however this ID integral only needs to be done once mak- 
ing this semi-analytic method for the image plane inte- 
gration computationally orders of magnitudes faster than 
classical ray shooting methods for computing the magni- 
fication of an extended source. 

4.3. Numerical Source Plane Integration 

Since the magnification must be computed numerically, 
the source plane integral in eq. |T|) must be computed nu- 
merically as well. The numerical integration scheme we 
use incorporates multiple grids with adaptive resolution 
to tile the source plane efficiently; details are given in the 
Appendix. 

5. RESULTS 

5.1. Three Probabilities 

As discussed in $2j we use a magnification cut to de- 
termine whether a system is labeled a strong lens or not. 
Since the labeling is applied only after follow-up observa- 
tions, we apply the cut to the total magnification (not the 
fiber magnification). We choose as our fiducial thresh- 
old ^ cut = 2, because this corresponds to the magnifi- 
cation of a point source located on the boundary of the 
multiply- imaged region for an isothermal sphere lens. A 
finite source magnified by this amount should show clear 
signs of strong lensing. In other words, we define 



to be the probability that a galaxy has a rogue emission 
line due to a source that is lensed. 

We should not just discard systems with magnifications 
below the cut. There is a range of positions (roughly 
speaking, behind the galaxy but outside the Einstein 
cone) where a source could send enough flux down the 
fiber to create a rogue line without being lensed. Such 
systems represent false positives in spectroscopic lens 
searches, and in order to understand the efficiency of 
a survey we need to assess the false positive rate. We 
define the complement of eq. ([23)1 to be the probability 
that a galaxy has a rogue emission line due to a source 
that is not lensed: 



P N (G;R src )=P G (G | 



(24) 



We can then let Pr = Pl + Pn be the total probabil- 
ity that a galaxy has a rogue line, while Rf — Pn/Pr 
is the false positive rate (defined to be the fraction of 
candidates for which follow-up observations do not show 
substantial evidence for lensing). Throughout our anal- 
ysis we keep track of the total rogue line probability Pr, 
the lensing probability Pl, and the false positive rate 
Rf. 

5.2. Back of the Envelope Estimate 

Before giving our full results, we can make a simple 
estimate of the lensing probability. For this estimate 
we ignore finite source and fiber effects and just take all 
sources inside the Einstein cone to be lensed. For simplic- 
ity we take the luminosity threshold Lq to be indepen- 
dent of source redshift, and we consider a non-evolving 
source luminosity function, so n s is constant. If we let Lq 
be the lowest luminosity threshold (i.e., computed from 
the flux limit So at the minimum source redshift), we 
should overestimate the number of detectable sources. 
With these simplifications, the lensing probability is 



Piest 
L ' 



n s Tr(D^b) 2 dD G 



D 



D 



3 Vc 
16tt 3 /ct\ 4 



□ V c max -A c ) 3 - 



16tt 3 /ct\4 



(!) 



P L (G;R S1 - C ) = P G (G | 

(^sre) ^* \^cut)- 



(23) 



° 3 \c 

[(1 + ^.max^.max - (1 + *i)A] 8 , (25) 

where C indicates a comoving distance. We take z SjIIlax = 
0.9. For fiducial LRG parameters zi = 0.2 and a = 250 
km/s, our estimate yields P£ st ~ 3.1%. In other words, 
we expect the lensing probability to be at the percent 
level. 

5.3. Dependence on Lens Parameters 

We first seek to understand how the various probabili- 
ties depend on the lens galaxy properties. Figured] shows 
the rogue line probability, lensing probability, and false 
positive rate as a function of source size, for different 
values of the lens redshift zi and velocity dispersion a. 
The total probability Pr for detecting a rogue line in an 
LRG spectrum is at the level of 2-3% (consistent with 
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FlG. 4. — Total rogue line probability (Pr, left), lensing probability (Pl-, middle), and false positive rate (Rp, right), as a function of 
the size of the source. We take as fiducial parameters (jgj, a, e, / y, 0») = (0.2,250 km/s, 0,0,0). In the top row we vary the lens redshift, 
in the middle row we vary the velocity dispersion, and in the bottom row we vary the velocity dispersion without including the effects of 
seeing. The spectroscopic fiber cut produces a non-monotonic dependence of Pr on a. This behavior also depends on the seeing. The 
lensing probability Pl turns over at large R STC (most visible in the a = 150 km/s curve) due to the magnification cut. The false positive 
rate also exhibits complex behavior with R S rc and is very sensitive to both a and Z\ . 



our estimate above), and depends moderately on lens 
redshift. Once a rogue line is detected, the probability 
that the line corresponds to a source that is significantly 
distorted decreases significantly with lens redshift: the 
false positive rate is only Rp ~ 35% at z\ — 0.15, com- 
pared with Rp ~ 80% at z\ = 0.45. The trend is not 
surprising because there is a finite upper limit on the 
source redshift, and the Einstein radius shrinks as the 
lens moves out in redshift (b oc D% 8 ). 
As expected, Figure [4] shows that the lensing probabil- 



ity depends strongly on a: Pl varies by almost a factor 
of 10 over the range 150 km/s < a < 400 km/s. For 
most values of a, the lensing probability increases with 
source size over the range < R SYC < 10 kpc. How- 
ever, for a — 150 km/s the curve reaches a peak at 
R src pa 5 kpc and then turns over. We attribute this to 
finite source effects: at a — 150 km/s the lens is simply 
not strong enough to significantly perturb a large source, 
(.Rein = 0.48 arcsec for z s — 0.9). In fact, all of the Pl 
curves would turn over if we went to large enough source 



size. 

We also find the surprising result that the rogue 
line probability does not increase monotonically with a. 
With 2 arcsec seeing, the ordering of the curves in in- 
creasing P R is a = 300, 250, 350, 200, 150, and 400 km/s 
(for all source sizes). In other words, the total rogue line 
probability for a = 150 km/s exceeds that for all other 
cases except a = 400 km/s. We attribute this to the fi- 
nite size of the spectroscopic fiber. When a is small, the 
Einstein cone is small and most sources inside the fiber 
are not lensed (indeed, the false positive rate is high) . As 
a increases, lensing can push some of the light outside 
the fiber (see Figure [2]), reducing the fiber flux and hence 
making the rogue line harder to detect. However, if the 
line is detected the chance that it corresponds to a lens 
is high. As a increases still further, the Einstein cone 
begins to fill the fiber (the false positive rate becomes 
very low), and the rogue line probability increases with 
a just like the lensing probability. 

The total rogue line probability, and its non-monotonic 
dependence on a, are sensitive to seeing. Eliminating 
seeing (bottom row of Figure [4]) changes the ordering of 
the Pr curves with different values of a; in particular, 
it shifts the case with the lowest rogue line probability 
to lower a. By contrast, seeing has less effect on the 
total lensing probability. These features can again be 
understood in terms of fiber effects. When a is low, the 
Einstein radius is small and most of the "lens" configura- 
tions lie well within the spectroscopic fiber, so seeing has 
little effect on the fiber flux and hence on Pl . Seeing can 
pull flux into the fiber from sources that lie outside, but 
these are predominantly non-lens configurations. The 
net result is that seeing increases the total rogue line 
probability, mainly by adding false positives. When a is 
high, by contrast, the fiber is mostly filled with "lens" 
configurations, so the false positive rate is low both with 
and without seeing. 

Though not shown here, we have also studied the ef- 
fects of varying ellipticity and shear parameters. We find 
no significant change in the lensing probability due to 
e and 7, in contras t to point source lensing statistics 
(|Huterer et al.l [20051 ) . The difference is presumably re- 
lated to the different statistical question (i.e., the proba- 
bility that a galaxy is a lens rather than the probability 
that a source is lensed), and to effects like the fiber cut 
that are specific to spectroscopic surveys. 

5.4. Dependence on Survey Parameters 

It is important to understand the dependence of the 
probabilities on our choice of the definition of a lens (the 
magnification cut /i cu t), and on the survey noise floor 
Sq. Figure [5] shows Pr, Pl, and Rf for 1.5 < /x cut < 3.0. 
Since /z cu t simply indicates whether a source is classi- 
fied as lensed or not, the rogue line probability does not 
change. However, the strong lensing probability Pl does 
depend on /x C ut- For i? src = 0.5 kpc, increasing /x cut 
from 1.5 to 3.0 decreases the lensing probability by a 
factor of 2.4. At R src = 10 kpc, the change is a factor of 
4.1. The dependence on source size can be understood 
in terms of the magnification regions in the source plane. 
As ^ cut is increased, the total area of the source plane 
with jj, > /x cu t decreases more slowly for small sources 
than for larger sources. 

If we double the noise floor Sq, we find that the 



rogue line probability is decreased by a factor of 0.65 for 
i? src = 0.5 kpc, and 0.63 for i? src = 10 kpc. This change 
is mainly caused by the number of sources brighter than 
the flux limit. Since n s oc T[l + a, Lo/(iL*\ and Lq oc 
So, doubling the noise floor changes the number of de- 
tectable sources by a factor of ~ r[l + a, 2Lo//xL*]/r[l + 
a, L / fiL*]. This works out to be a factor of 0.6-0.8 for 
typical Lq and fi values. This simple estimate agrees 
quite well with our full calculations despite ignoring com- 
plicated seeing effects, magnification effects, and the red- 
shift dependence of Lq. 

5.5. Total Probabilities and Higher Order Statistics 

We are now ready to compute overall probabilities by 
averaging over an appropriate sample of galaxies (see eq. 
2]) . First we need to consider how many galaxies we need 
to include to obtain accurate statistics. Figure E] shows 
the three probabilities as a function of the number of 
LRG realizations (for R S1C — 0.5 kpc and (3 = 0, meaning 
no evolution in the source [O II] luminosity function). 
All three probabilities converge rather quickly, indicating 
that our fiducial sample size of 800 LRGs is sufficient. 

Figure [7] shows the total probability for detecting rogue 
emission lines in LRG spectra, the probability of identi- 
fying lensing in follow-up observations, and the SLACS 
false-positive rate all as a function of source size averaged 
over our full sample of 800 galaxies. 

At this point we also examine the effects of evolution 
of the source [O II] luminosity function. We consider 
four models, ranging from no evolution (/3 = 0) to a 
strong increase in the number of [O II] emitters in the 
past (0 = 3). A value of sa 2 -3 for z < 1 is pre- 
ferred by iGlazebrook et all (|2003f ) based on SDSS de- 
rived star formation histories. LF evolution tends to 
decrease the rogue line and lensing probabilities. The 
explanation is that we fix the total number density of 
sources in the range 0.35 < z < 1.5 from the observa- 
tions bv lHogg et al.l (|1998l ). so increasing shifts a higher 
fraction of those sources to higher redshifts, and there- 
fore decreases the number of sources below the SLACS 
upper limit z S!inax = 0.9. 

For a source size of 0.5 kpc, we find the total probabil- 
ity for rogue line detection varies from 2.0% for source 
LF's with no evolution {0 = 0) to 0.9% for LF's with 
strong evolution [0 = 3). With i? src = 10 kpc, the prob- 
abilities are 3.0% (f3 = 0) a nd 1.3% (0 = 3). Th e impli- 
cation is that in the original [Bolton et al.l (|2004[ ) sample 
of 50,996 LRG spectra, ^460-1,530 should contain rogue 
[O II] emission lines in their spectra above a threshold of 
6 x 10 -17 erg/s/cm 2 . Of those, ~250-640 should show 
significant strong lensing features when imaged with high 
spatial resolution. The remaining systems are false posi- 
tive detections in which there is a background galaxy that 
is not significantly perturbed by lensing effects. The to- 
tal false positive rate Rf{R S ic) also varies significantly 
with with source size and but has typical values Rf ~ 
50% (see Figure [JJ . It is important to note that while 
the breakdown into lenses and false positives depends on 
our choice of the magnification cut (/i C ut = 2.0), the total 
number of rogue lines is a robust prediction. 

While our fiducial results have been computed for 2 
arcsec seeing, it is instructive to consider how seeing af- 
fects our result. Figured] shows the three probabilities as 
a function of source size for various values of the seeing. 
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The total number of rogue lines remains unchanged (of course), but the lensing 
probability Pj, varies by a factor 2.4 at small source size and 4.1 at large source size. As expected, the false positive rate Rp is sensitively 
dependent on the definition of a lens, /t C ut- 
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Fig. 6. — Total detection probabilities as a function of the 
number of LRG realizations. In this convergence test we fix the 
source size to be R BTC = 0.5 kpc and we use a non-evolving source 
luminosity function (J3 = 0). The probabilities converge quickly, 
indicating that 800 realizations is more than adequate to yield ac- 
curate statistics. 

The effect of seeing is to introduce a "tilt" to the lens- 
ing probability curve, giving a moderate boost at small 
source size and reduction at large source size. Seeing 
increases the total rogue line probability, especially at 
small source size. The implication is that seeing en- 
hances the number of detected rogue lines that do not 
correspond to lens systems at all source sizes but most 
dramatically at small R src . Indeed, the false positive rate 
increases monotonically with R src when seeing is unim- 
portant (the FHWM is small compared with the size of 
the spectroscopic fiber), while it flattens out and devel- 
ops a minimum near R src 5.5 kpc when the seeing is 
important. 

Although there is significant uncertainty due to source 
size and LF evolution effects, our lensing estimates are 
clearly higher than the ~150 rogue lines and ~60 lens 



systems found in the initial SLACS sample. The most 
likely explanation for this discrepancy is that the SLACS 
selection criteria require two additional emission features 
(HP, [O III] 4959, or [O III] 5007) besides the primary 
[O II] 3727 line. The presence of these additional fea- 
tures was required to substantially reduce the number 
of false positives, but may have eliminated many real 
lens systems as well. Incorporating multiple emission 
line statistics into our calculations would require knowl- 
edge of the joint probability distribution between [O II] 
and secondary line luminosities and is beyond the scope 
of the present study. Furthermore, secondary lines tend 
to appear in the ~ 7000-9000 A region of the spectrum 
where sky noise is more problematic; many spectra with 
secondary features buried in the noise may have been re- 
jected as targets. In any event, we predict that there is 
a large number of g-g lenses waiting to be discovered in 
the SDSS spectroscopic data. 

Finally, it is interesting to consider how the galaxies 
that produce g-g lenses may form a biased subset of all 
LRGs. Figure [H] shows the distributions of redshift and 
velocity dispersion for our parent population of galaxies, 
and for the same population with each galaxy weighted 
by its lensing probability (with R SIC — 0.5 kpc, p = 
0, seeing=2"). There is clearly a bias towards larger 
velocity dispersions, which is natural, but it is weaker 
than in traditional quasar lens statistics because of fiber 
cut effects. There is also a bias towards lower redshifts, 
because lower-redshift galaxies have larger Einstein cones 
(cf. Figured]). We test the null hypothesis that the two 
populations represent the same prob ability distr i butio n 
with the Kolmogorav-Smirnov test (|Press et alj|1992j) . 
Comparing the distributions of the velocity dispersion 
(redshift) yields only a 2.5 x 10~ 7 % (0.05%) chance that 
the populations are drawn from the same distribution. 

6. CONCLUSIONS 

We have introduced a statistical method to estimate 
the expected number of strong galaxy-galaxy lens sys- 
tems in a spectroscopic survey. In the process, we have 
also developed a semi-analytic technique for calculating 
the magnification of a finite source due to an isothermal 
ellipsoid galaxy in an external shear field. Conceptually, 
the most important result is that the size of the spec- 
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Fig. 7. — Total detection probabilities after averaging over 800 sets of LRG parameters, for four different models of LF evolution. Given 
these results, we expect that the initial SLACS sample of 50,996 LRGs should contain ~460-l,530 galaxies with rogue emission lines in 
their spectra, and ~250-640 should reveal strong lensing features in follow-up observations. 
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Fig. 8. — Total detection probabilities as a function of source size and the FWHM of the Gaussian seeing model ({3 = 0). Seeing effects 
boost the lensing probability at small -Rare but suppress the lensing probability at large -R src . The total rogue line detection probability is 
boosted for all source sizes. The fraction of non-lenses is much more sensitive to seeing effects at small -Rsrc leading to the minimum in the 
false positive rate in Figure [7] 



troscopic fiber provides a significant selection effect. In 
our analysis of the SLACS survey, the fiber cut yields 
the unexpected result that the probability of detecting 
a rogue emission line is a non-monotonic function of the 
velocity dispersion a. Since larger a corresponds to a 
larger Einstein cone, one would naively expect the rogue 
line probability to increase with a. However, large a also 
corresponds to a large image separation, which can cause 
much of the source flux to fall outside the spectroscopic 
fiber. The situation is further complicated by the effects 
of atmospheric seeing which can add flux into the fiber 
from images outside or remove flux from the fiber from 
images inside. It will be crucial to compute the effects of 
the fiber cut and seeing, customized to the parameters 
of the spectrograph, for any future search for g-g lenses 
in large spectroscopic surveys. 

The lensing probability has a fairly weak dependence 
on the size of the source, but a stronger dependence on 
the evolution of the source luminosity function. Lensing 
introduces biases such that lens galaxies will tend to have 
higher velocity dispersions and lower redshifts compared 
with the parent population of galaxies. 



Incorporating parameters from the SLACS survey, we 
estimate that there should be ~460-l,380 rogue emission 
lines in the sample of 50,996 LRG spectra. Of these, 
~250-640 should show clear evidence of strong lensing 
in follow-up observations. The broad range of proba- 
bilities is due primarily to uncertainties in the physical 
size of [O II] emission regions and in the evolution of 
the [O II] luminosity function. Specifically, small sources 
and strong evolution yield lower probabilities while large 
sources and no evolution give higher probabilities. Our 
estimates are notably higher than the numbers actually 
observed in the SLACS survey so far. We attribute this 
to their requirement that multiple emission lines be de- 
tected, to the large amount of sky noise in the long wave- 
length region of the spectra, and to the potential extinc- 
tion of the [O II] emission line by dust in the lens galax- 
ies. Our calculations imply that there are many more 
galaxy-galaxy strong lenses waiting to be found in the 
SDSS spectroscopic data. 

While our methods have been specifically applied to 
the SLACS survey, they should be applicable to all finite 
source lens searches in upcoming spectroscopic surveys. 
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Fig. 9. — Histograms of velocity dispersion (left) and redshift (right), for the LRG parent sample (solid line) and the subsample that 
represent lenses (dashed line, R SIC = 0.5 kpc, /3 = 0, seeing=2"). All histograms are normalized to unit area. The lens galaxies are biased 
towards higher velocity dispersion and lower redshift. A KS test comparing the distributions of the velocity dispersion (lens redshift) yields 
only a 2.5 X 10 — 7 % (0.05%) chance that the populations are drawn from the same distribution. 



Future improvements to our method would involve in- 
corporating non-uniform brightness distributions for the 
background source galaxies. 
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APPENDIX 

TWO-DIMENSIONAL INTEGRATION 

Our source plane integration is carried out via a 2-D integrator that tiles the integration region efficiently by using 
a movable grid with adaptive resolution. As an example to illustrate our method, we first consider the magnification 
calculation for a uniform brightness finite source (with size R src centered at the origin) which is lensed by an isothermal 
ellipsoidal mass distribution (SIE). 



12 





Fig. 10. — A schematic representation of the tiling algorithm. The algorithm uses a macrogrid with adaptive resolution that subdivides 
if the number of corners that are "in" is 7^ or 4. Once the maximum number of sublevels is reached, the macrogrid is shifted in a spiraling 
geometry about the origin. 



Since the magnification is denned as, 

Area of the images , . , 

= -A r — i — , ( A1 ) 

Area of the source 

we must find the total area of the images by integrating over the 2-D image plane. We start with a single "macro" -grid 
whose lower right corner is centered on the origin and suppose that the edge of an image passes through this grid as 
shown in Figure [TO] A grid point (x,y) is flagged as "in" if, 

[x - a x (x,y)] 2 + [y- a y (x,y)} 2 < R 2 rc (A2) 

and "out" if, 

[x - a x (x,y)} 2 + [y- a y (x,y)} 2 > R 2 IC (A3) 

Here a(x, y) is the deflection angle due to the SIE at point (x, y) and relates the image plane grid points to the source 
plane coordinates (u, v) via the lens equation, 

u = x - a x (x, y) 

v=y - u y {x,y) (A4) 
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Fig. 11. — A schematic representation of the evaluation order for the "steps" for the spiraling algorithm. The origin is marked with an 
O and the spiraling terminates once N to [ revolutions yield no significant change in the integrated area. 



We then recursively subdivide portions of this grid until all four corners of each subgrid are either "in" or "out" (or 
until the maximum number of subdivisions has been reached). For a given subgrid this implies that our criterion for 
further subdivision is: 

Ni n = 1, 2, or 3 : subdivide 

Ni n = or 4 : do not subdivide ^ ' 

where Ni n is the number of corners that are flagged as being "in" the image. 

Once this initial grid has been appropriately subdivided, we then move the macrogrid in a spiraling geometry and 
repeat the process (see Figure 111 [) keeping a running total of the image area upon each revolution of the spiral. This 
spiraling procedure self-terminates once the macrogrid has completed N to i revolutions with little to no change in the 
total image area. 

Our technique has the advantage of achieving high resolution with relatively few subdivisions per grid and, more 
importantly, does not require preset boundary conditions. That is, the self-termination criterion simply stops the 
spiraling procedure once the total image area is no longer changing. However, care must be taken to avoid missing 
"distant" images (images which appear far from the origin) , and in practice we do set a minimum number of revolutions 
for the spiral. 

In the previous example, we assumed our function evaluations at the grid points took only two values: 1 or (i.e., 
"in" or "out" of the image). The next step is to allow each grid point to assume a continuous value f(x,y). With this 
generalization, we must modify our subdivision criterion which we now take to be, 

\F avg -F c \ > e\F avg \ : subdivide 

\F aV g—F c \ < e \F avg \ : do not subdivide, ,; 



where 



Favg = ^[/(xi,yi) + f(x 2 ,yi) + f(xi,y 2 ) + f(x 2 ,y2)] (A7) 



4 

is the average of the function values at the four grid corners and 



F c = f[ X -^^\ (A8) 



is the function value at the center of the grid. The tolerance value e must be set by hand and is roughly a measure of 
how much f(x,y) varies over the region [x\ : atra] [z/i : V2\- That is, our grid only increases resolution in regions where 
the function varies rapidly with either x or y. Since the spiral termination criterion is still such that the total integral 
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does not significantly change after N to i revolutions, this 2-D integration scheme is ideally suited to integrate functions 
for which f(x, y) — > as x, y — > oo. 



